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Electroweak transition form factors of heavy meson decays are important ingredients in the extraction of the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix elements from experimental data. In this work, within a light-front 
framework, we calculate electroweak transition form factor for the semileptonic decay of D mesons into a pion or 
a kaon. The model results underestimate in both cases the new data of CLEO for the larger momentum transfers 
accessible in the experiment. We discuss possible reasons for that in order to improve the model. 
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1. Introduction 

Charmed mesons inaugurated a new era in par- 
ticle physics. The discovery of the first heavy 
meson was an evidence in favor of the standard 
model as giving the correct description of low en- 
ergy particle physics. The new quark was con- 
firmed experimentally (see [1]) in a meson quark- 
antiquark state, i.e., in a heavy- light quark sys- 
tem. Since then, the study of semileptonic decays 
gained importance, mainly in the last years, as 
these decays constitute a laboratory to test the 
standard model predictions and to get informa- 
tion about the CKM [2] matrix elements. These 
matrix elements, Vcq, can be extracted from the 
experimental data of D decay. However, theo- 
retical inputs are necessary in order to interpret 
the data from the decay experiments as obtained 
by BES [3], Babar g]. Belle [5], and CLEO [617] 
collaborations. 

The QCD nonperturbative effects in the 
semileptonic decay are modelled with transitions 
form factors that have been calculated with lat- 
tice QCD jg, QCD sum rules [9] and light-front 
constituent quark models |10I11I12] . This mo- 
tivates us to study the semileptonic decay of 
D mesons within a light-front covariant model 
(LFCM) proposed in rcf. |13] and applied to cal- 



culate the pion, kaon and D electromagnetic form 
factors. 

In general, light-front models of hadrons make 
possible the calculation of many obscrvablcs, with 
the advantage of working at the amplitude level 
obtained directly from matrix elements taken be- 
tween hadron light-front wave functions. The va- 
lence component of the wave function allows a di- 
rect interpretation of the relevant model parame- 
ters, being also an useful concept when analytical 
forms of the hadron Bethe-Salpeter amplitude are 
the starting point to calculate the physical quan- 
tities. Applications of light-front models have 
been performed in many cases like in the calcula- 
tion of pscudoscalar meson decay constants |14j , 
electromagnetic form factors of the p-mcson |15j , 
pion jl6|17| and kaon [18]. In particular con- 
stituent light-front quark models (CLFQM), are 
successfully in describing the pion electromag- 
netic form factor compared with the new exper- 
imental data [19] (see e.g. [20] and references 
therein) . 

One important theoretical aspect for the cal- 
culation of electroweak form factors performed 
within the light-front approach, is the inclusion 
of the non-valence contribution besides the va- 
lence one when the Drell-Yan condition is not 
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respected, as in the case of the computation of 
time-hke from elastic and transition factors (see 
a discussion in 20 ). The nonvalence contribu- 
tion to the electroweak matrix elements of the 
current are associated with Fock-components of 
the wave function beyond the valence, that can 
be as well attributed to two-body currents acting 
on the valence wave function [ 21^^2] . Notably, 
in some particular light-front models, the non- 
valence contribution should be included in the 
calculation of electromagnetic current matrix el- 
ements to satisfy the requirements of covariance 
in the Drell-Yan frame TW. Moreover, in calcula- 
tions of electroweak observables within light-front 
frameworks, apart from the frame used, it is im- 
portant to note that the contribution of nonva- 
lence terms depends on which component of cur- 
rent is employed to extract the form factors (see 
[I6|23] for details). 

Recently, we have proposed an analytical co- 
variant model of the Bethe-Salpeter amplitude 
of mesons [13] and we have used a Wick rota- 
tion of the light-front energy in the Mandelstam 
formula to calculate electromagnetic form fac- 
tors (LFCM). The Wick rotation applied to the 
minus component of the loop-momentum allows 
a direct four-dimensional integration avoiding the 
singularities of the Bethe-Salpeter amplitudes in 
Minkowski space in the calculation of the observ- 
ables. The computation takes care implicitly, all 
possible zero-modes or nonvalence contributions 
to the matrix elements of the electromagnetic cur- 
rent. We have used the above model and calcu- 
lational strategy to obtain electromagnetic form 
factors in the space-like regime and also decay 
constants of light and heavy pseudoscalar mesons. 

Once elastic form factors of pseudoscalars have 
been calculated in the model, it is worthwhile to 
pursue the investigation of transition form factors 
in heavy meson decay. The semileptonic decays 
of the heavy mesons offer an useful laboratory to 
test theoretical models of the corresponding tran- 
sition form factors in the weak sector of the stan- 
dard model. Here, the covariant Bethe-Salpeter 
vertex model and computational strategy using 
Wick-rotation on the light-front are applied to 
calculate the electroweak transition form factors, 
f+{q^) and /_ (q^) for D meson decays into a pion 



or a kaon. In sect. 2, we describe the model and 
give their parameters. In sect. 3, the form factors 
and kinematics for the decay are presented. The 
numerical results and conclusions are presented 
in sect. 4. 

2. The vertex model 

The vertex model of the Bethe-Salpeter ampli- 
tude for the meson-g— q proposed in reference [13] 
is given by the following expression: 

AM{k,p) = {e-ml){{p-kf-ml) (1) 

Tm 

""(fc^-Ai.+^e)" ((p-fc)2-A2, + ze)"- 

Here, Am is the scale associated with the me- 
son light-front valence wave function and n is the 
power of the regulator. The masses of the con- 
stituent quarks within the meson state are uiq 
and mg. In the vertex function, the factors in 
the numerator avoid the cut due qq scattering 
if ruq + rnq < Mh- In order to confine the 
quarks, the regulator scale \m is chosen such that 
\m > Mh/2, where Mh is the meson mass. This 
constraint does not allow the meson decay into 
particles of mass Xm ■ 

The general matrix elements for the electro- 
magnetic current, are written as a tree-point func- 
tion, using the Mandelstam formula, given in the 
equation below: 

(p' \Ji; {q^)\p)^-^[d'k Tr[0]F{k,p,p').i2) 
[2tt) j 

where the operator O is 

O - [Sf {k - p') rSpik - p)Sf (k)] . 

NcSpip) is the Feynman propagator for con- 
stituent quarks and 

F{k,p,p') = Am' {k,p')AM {k,p). 

The electromagnetic form factor for pseu- 
doscalar mesons is obtained from the matrix ele- 
ment of the conserved vector current as: 

{p'\j!;{q')\p) = {p + prF^^Si<i'), (3) 

with the choice of plus component of the current, 
J+ = + J^, that has the minimal contribu- 
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0.345 
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0.346 


2138 


220 
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Table 1 



LFCM parameters for the pion, kaon and . 
For each meson the range parameter Am of the 
Bethe-Salpeter vertex ^ with rt = 10 is fitted 
to the corresponding pseudoscalar decay constant 
for the given constituent quark masses. The val- 
ues used are = 94 MeV, Jk = 113 MeV and 
fo = 144.5 MeV [I]. The different constituent 
quark masses originate the parameterizations la- 
belled by P, K and D, for the pion, kaon and D, 
respectively. 

tion from the instantaneous terms of the Feyn- 
man propagator of the constituent quark. Gener- 
ally, this "good" component of the current is pro- 
tected against the contribution from zero-modes 
in the Drell-Yan frame, although a careful anal- 
ysis has to be performed for each choice of the 
vertex, when the integration over the minus com- 
ponent of the loop-momentum is done analyti- 
cally |16|23j . The normalization condition of the 
vertex function is such that i^p™(0) — 1. 

The weak decay constant of the pion is defined 
by the matrix element of the axial-current: 

(0|A^(0)b) =^V2/psP^ (4) 

where = g(x)7^7^^g(a;), with analogous ex- 
pression for the kaon and D. 

The final expression for decay constant of the 
pseudoscalar mesons is given by 

^TrmAM{k,p), (5) 
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Figure 1. Transition form factor f+{q^) normal- 
ized to /+(0) = 1 for Z) — >■ ttIv. Calculations with 
the LFCM model compared to results from lattice 
[25] and dispersion relation calculations |26) . 

where PS ^ (7r,K,D) and 

g = Tr[7V5f(fc)7'"^5F(fc-p)] . 

The plus component of the axial current is 
adopted as discussed above for the electromag- 
netic current. The numerical integration in the 
Euclidean minus moment component allows one 
to use, equally well, any other component. 

The vertex model parameters for the pion, kaon 
and D are obtained from the decay constants for 
given constituent quark masses as presented in 
table 1. Our choice of constituent quark masses 
follows the suggestion that the chiral symme- 
try breaking mechanism produces the light con- 
stituent quark mass, and the other quarks like s 
and c have in addition the current quark mass. 
In our case we use the current quark masses of 
124 MeV and 1239 MeV for s and c, respectively. 
We follow the reasonings presented in [53], and 
the current quark masses are within the accepted 
range of values (see [T]). 

3. Transition Form Factors 

The transition form factors of the weak vector 
current investigated here appears in the semilep- 
tonic processes D — nlv and Ds — > Klv de- 
cays. In particular, the leptonic decay of Dg 
allows to extract the decay constant. In the 
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Figure 2. Transition form factor f+{q^) nor- 
malized to /+(0) = 1 for I? — > Klv. Calcula- 
tions with the LFCM model compared to results 
from lattice and dispersion relation calcula- 
tions pS] , 

semileptonic decay of a heavy pseudoscalar me- 
son to a light one, the nonperturbative dynamics 
of QCD appears through the matrix elements of 
the weak vector current parameterized in terms of 
two transitions form factors, f+{q^) and f-{q^): 

{7:{K)\V'^\D{D,)) = iP' + Prf+{q^) + f^{q^)q^, 

where = qji_iC, with q — (d, s), and = 
(P' — P)'^. Pf^ stands for the four- momentum of 
the D meson and P'^ for the light pseudoscalar in 
the final state. The matrix element of the weak 
vector current can be written also as: 



{7r{K)\V^\D{D,)) 



{P' + P) 



'foil 



(6) 



when the momentum transfer vanishes /+(0) ~ 
/o(0). In the calculation of transition form fac- 
tors, we use the plus component of the weak vec- 
tor current. 

The reference frame used in the case of the 
elastic electromagnetic form factors is defined by 
[US]: P+ = P'+ + q+, <z+ > 0, Pi = Pi = 



and P' = P^ = M^. For the decay process, we 
use the following reference frame: P^ = P'^—q^, 
q^ < and P^ = P± = 0. The numerical inte- 
gration of the matrix elements is done with the 
Wick rotation of the minus momentum as formu- 
lated in 113]. 

4. Results and Conclusion 

In the LFCM models used use n = 10, we have 
chosen the range parameters A m by fitting the ex- 
perimental pseudoscalar decay constants, the are 
presented in the table 1, together with the charge 
radius. The experimental pion charge radius is 
0.672 fm[Tn] much larger than our results, that 
comes as a consequence of the regulator power 
n = 10. The same is observed for the kaon which 
has an experimental charge radius of 0.56 fm |19) . 
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BES [3, 


0.780± 


0.730± 
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0.695±0.023 


0.624±0.036 


BaBar [5] 


0.728± 




CLEO [6] 


0.739±0.009 


0.665±0.019 


CLEO-c [2] 


0.763± 


0.628± 


LQCD [g 


0.733± 


0.643± 


Model 






(K1,D1,-) 


0.729 




(K2,D2, - ) 


0.772 




(K3,D3,P3) 


0.913 


0.491 


(K4,D4,P4) 


1.020 


0.739 


(K5,D5,P5) 


1.100 


0.936 



Table 2 

Experimental, lattice and model transition form 
factor results for q^ — 0. 

We should observe that smaller power in the 
vertex function AA/(fc,p) give a good fit for low 
energy observables calculated with this model [13] 
and describe well the electromagnetic form fac- 
tors of the pion and kaon jl3i . But, in the case 
of the transition electromagnetic form factors we 
need the higher power n, in order to reproduce 
the experimental transition form factor /^(g^) in 
the low momentum region (see figures 1 and 2). 
While the slope near q^ = seems reasonable, 
the calculation underestimate the experimental 
data for large momentum transfers for the range 
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of constituent quark masses used. The effect of a 
light quark mass variation of 10% is seen mainly 
at large momentum transfers. 

In table [51 we show transition form factors 
/+(0) for D and Ds compared to the experimen- 
tal data 1614151713] and to lattice calculation [25] . 
The sensitivity to the constituent quark mass is 
strong, and no simultaneous fit of D and Dg tran- 
sition form factor at zero momentum transfers 
was achieved by the model. The small sizes of the 
pion and kaon are presumably responsible for the 
missing curvature demanded at large momentum 
transfers by the data. The large power implied in 
a high value of the scale Aj\/ of about 3 GeV for 
the light mesons, to compensate for the power 
n = 10 while fitting the decay constant. This 
hard scale as seen in table 1 shrinks both pion 
and kaon, as well as the D meson. The overlap 
between the wave function represented by /'*'(0) 
and the slope of the transition form factor seems 
not that sensitive to these hard scales but to the 
physical decay constants. Therefore, the vertex 
has to have a much softer dependence in momen- 
tum, and a smaller power should be required for a 
better fitting of these observables. Indeed the ex- 
pected behavior of the vertex function with mo- 
mentum should have a power about n = 2, as 
expected from the one-gluon exchange and the 
monopole decay of the pion electromagnetic form 
factor. Therefore, it is expected that even the 
D meson should be somewhat larger than the 
present calculations indicate of a charge radius 
of 0.34 fm, in order to accommodate a reasonable 
fit to the transition form factors. 
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